Thermally-induced vacuum instability in a single plane wave 
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Ever since Schwinger published his influential paper [J. Schwinger, Phys. Rev. 82, 664 (1951)], it 
has been unanimously accepted that the vacuum is stable in the presence of an electromagnetic plane 
wave. However, we advance an analysis that indicates this statement is not rigorously valid in a real 
situation, where thermal effects are present. We show that the thermal vacuum, in the presence of 
a single plane- wave field, even in the limit of zero frequency (a constant crossed field), decays into 
electron-positron pairs. Interestingly, the pair-production rate is found to depend nonperturbatively 
on both the amplitude of the constant crossed field and on the temperature. 

PACS numbers: 12.20.Ds, ll.10.Wx 



It has long been known that an inevitable conse- 
quence of Dirac's theory of the electron is that in re- 
gions of sufficiently high energy density, the quantum 
vacuum can break down in a spontaneous generation of 
electron-positron pairs. Following the initial results of 
Sauter Heisenberg and Euler Q and Wcisskopf 
in a seminal work, Schwinger Q derived a central re- 
sult of strong-field quantum electrodynamics, the rate 
per unit volume of pair creation R in a constant and 
uniform electric field of strength E, of leading order be- 
haviour R = (E/E CI ) 2 (c/X 4 )(Sir 3 )- 1 exp(-TrE cr /E), for 
E/E CI <C 1, positron charge e, mass m, Compton wave- 
length X = h/mc and so-called "critical" electric field 
E cr = m 2 c 3 /eh = 1.3 x 10 16 Vcm" 1 . A combination 
of factors has intensified research efforts to better under- 
stand and devise ways of detecting this phenomenon. On 
the one hand, upcoming laser facilities areplanncd, such 
as ELI (Extreme Light Infrastructure) [5| and HiPER 
(High Power laser Energy Research) [6(, that intend to 
reach fields as large as a percent of the critical value. On 
the other hand, recent theoretical results, both analytical 
[7 - 11 1 and from numerical simulation [lij , strongly indi- 
cate that at orders of magnitude well below the critical 
field, pair creation could be observed jl~3j ] . 

Decay of the vacuum is predicted to occur in a vari- 
ety of contexts. In intense electromagnetic backgrounds 
such as two plane waves propagating in different direc- 
tions lL4H16| or a plane- wave and Coulomb field combina- 
tion [l7| , but also in more exotic contexts, such as being 
seeded by magnetic fields of magnetars |l8( . Another ex- 
ample is through thermal radiation, with pair-creation 
rates having been calculated in a constant electric field 
in various formalisms l3"21|, also including stimulated 



pair creation [22} , However, it has long been accepted 
that pair creation can never occur in single plane waves, 
as encapsulated in Schwinger's famous statement "there 
arc no nonlinear vacuum phenomena for a single plane 
wave, of arbitrary strength and spectral composition" 
3- The physical reason for this statement is that all 
photons in a plane wave propagate in the same direction 
and so cannot interact with each other. The mathemat- 
ical origin lies in the fact that the two electromagnetic 
invariants vanish for a single plane wave [i| . 

In the current letter, we seek to demonstrate how pair 
creation can, in fact, proceed, when the vacuum is po- 
larised by a single plane wave, if one acknowledges the 
inevitable presence in all real physical scenarios, of back- 
ground heat radiation. We show that this is the case even 
if the frequency of the plane wave tends to zero (the so- 
called "constant crossed field" configuration). Moreover, 
by considering the existence of a thermal background, we 
will derive an expression for the rate of real pair creation 
which is non-perturbative in both the field strength E 
and the background temperature T, in both prefactor 
and exponent. Apart from the conceptual advance, our 
results may be of relevance for typical applications of the 
Schwinger effect such as string-breaking models [23| for 
heavy-ion collisions or for a full understanding of the pair 
production processes near neutron stars. In those and 
other cases, the thermal bath is an important ingredient. 

We consider a scenario in which a region of the vac- 
uum polarised by an external electromagnetic field can 
be considered as bathing in a thermal background, repre- 
sented by a photon gas in equilibrium, at temperature T. 
The high external field intensities that interest us can be 
characterised with the classical non-linearity parameter 
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FIG. 1. The envisaged scenario of thermal photons (wavy 
lines) interacting via the polarisation operator (thick circles) 
in equilibrium with an intense external constant crossed field 
(E, B) (thicker gray lines, indicating the direction n = E x 
B/EB). 



£ = (m I 'uji)(E I ' E a ) 3> 1, with uji the angular frequency 
of the external field photons, where we set here and sub- 
sequently h = c = kB = 1- For £ 3> 1, pair-creation from 
vacuum proceeds mainly via tunneling and in such pro- 
cesses, the frequency of the external field ceases to play 
a role, with the rate tending to that in which the limit 
lui — >• is taken [24[. For a plane wave, this describes 
the "constant crossed field" background, where the elec- 
tric field E and the magnetic field B are constant, equal 
in strength (E = B) and perpendicular to one another. 
Noting that modern laser systems already allow for val- 



ues of £ of the order of 10 [25j, we choose to study the 



behaviour of a photon gas in thermal equilibrium with a 
constant crossed field. 

The leading-order interaction between thermal pho- 
tons and the external field is contained in the polari- 
sation operator, as displayed in Fig. [T] In general, the 
polarisation of the vacuum in a background field F^(x), 
with four-co-ordinate x, which varies slowly over a Comp- 
ton wavelength X = 1/m, can be investigated via the ef- 
fective Lagrangian of QED (26j. It is this object that 
Schwinger referred to as being identically zero for all 
plane-wave backgrounds. In his analysis of a constant 
field and for the present case F^ v {x) = F^, there are 
only two relativistic invariants related to the background 
electromagnetic field F^: -F^F^ = 2(E 2 - B 2 ) and 
- F ^F* v = 4(E ■ B), where F* v = e^F^/2 with 
Cftvaf) being a tensor antisymmetric in all indices and 
E and B the electric and magnetic field. In a plane 
wave, and in particular in a constant crossed field, both 
the above invariants vanish identically. However, in the 
presence of a thermal bath, the virtual electron-positron 
pairs effectively couple the polarising background field 
with the thermal photons. If = (uj, k) indicates the 
on-shell four-momentum of a thermal bath photon, there 
exists a third relativistic invariant (k^F^) 2 , which is in 
general non-zero in a plane-wave background. This pa- 
rameter is related to the so-called quantum non-linearity 
parameter \ = e-\/|(fc /J i 7 ' AIl/ ) 2 |/m 3 . For photons scatter- 
ing in a constant crossed field (E, B) of unit wavevector 
h = E x B/EB, we have x = X-e( w — k ■ Je)/m, where 
Xe = E/E CT . Therefore, this will be the defining mi- 
croscopic variable describing the interaction between the 



photon gas and the external constant crossed field. After 
having averaged over the photon vectors fc M in the ther- 
mal bath, the additional relativistic invariant may also be 
written as (M^f") 2 , where u M is the 4-velocity vector of 
the heat bath. In the heat-bath rest frame, this invariant 
corresponds to — E 2 . In the following analysis, we limit 
ourselves to more accessible scenarios in which \ can be 
regarded as much less than 1. Qualitatively this occurs 
if the parameter 5 = (T/m)xE = (T/m)(E/ E cx ) is much 
smaller than 1. In the following, we provide additional 
conditions for the validity of our approach. 

In regarding the thermal vacuum polarised by the ex- 
ternal field as filled with a gas of photons, we are im- 
plicitly making two assumptions: i) that the polarisation 
from thermal photons can be included perturbatively and 
ii) the rate of pair creation from purely thermal effects 
is negligible. The first of these assumptions is valid if, 
in the formation volume of the pair production process, 
the number of thermal photons is less than unity. The 
formation length Lf, in which on avera ge a single pair is 
generated, can be seen to be Lf = X/xe [H|, whereas the 
density of photons in the gas is given by integrating the 
Bose-Einstein distribution over d 3 k to give 2£(3)T 3 /ir 2 , 
with £(•) the Riemann Zeta function. Therefore, the av- 
erage number of thermal photons in the formation vol- 
ume L 3 is 2(({3)/ir 2 ){x E T/m) 3 = 2(C(3) A 2 )<5 3 . So by 
fulfilling the condition S 3 -C 1, the first assumption will 
hold and in this limit, we can also obtain the total pair 
creation rate as an incoherent sum over the thermal en- 
semble. The parameters for which the second assumption 
is valid, can be found by first considering the case of a 
vacuum solely populated by a thermal photon gas. The 
leading order pair creation process then comes from a 
collision of two free photons [271 ] . which can be summed 
over the thermal ensemble to give the background pair 
density p2~ ( (T) = dN2 1 /dVdt for number of pairs cre- 
ated A^-y with volume V and time t. For T/m -C 1 this 
becomes: 
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where a is the fine-structure constant. This background 
rate will later be compared to the rate of pair produc- 
tion arising from the interaction of the thermal bath with 
the constant crossed field to show where also the second 
assumption, and hence our analysis, is valid. 

Bearing these remarks in mind, we can proceed via the 
optical theorem. This allows one to relate the imaginary 
part of the polarisation operator in an external electro- 
magnetic field to the rate of electron-positron pair cre- 
ation for a single photon propagating in this background, 
Ra{k,XE) = lm[l%n. lhV (k,XE)l* a v ]/u, for normalised pho- 
ton four-polarisation 1% (Y,Ja l T ~> -9^), a = {1,2} 
and polarisation operator Yi^ v {k,XE)- For brevity, we 
insert the derived polarisation operator in a constant 
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crossed field background (see e.g. 28[) into this expres- 
sion, giving a single-photon rate: 

2am 2 /"» 2« + l + 3(-ir Ai'(z) 

Ra{k,XE) = 5 / Cfo , --. , 2) 

where z = (v/x) 2 ^ 3 , Ai(-) is the Airy function of the first 
kind [29|. 

In the situation of interest where the value of the pa- 
rameter S is very small, we can regard the condition 
X <C 1 to be satisfied by all save a negligible number of 
photons. Therefore, we can asymptotically expand Eq. 
([5]), which will be further justified in the discussion, to 
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The density of pairs created by the external field-thermal 
bath interaction pth(T, xe) = dNth/dVdt, can then be 
obtained by summing Eq. ([3]) over the Bose-Einstein 
distribution, well approximated by an integration: 



Pth(T, X E) = Y. J (2^)3 ^Jt^ R ^Xe) 



(4) 
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FIG. 2. Logarithm of the expected number of pairs generated 
in a constant crossed field in the tiny space-time volume A 4 = 
7.4 x 10 -53 cm 3 s. The area above the solid line signifies where 
the background process of pair creation due to purely thermal 
effects (see Eq. {TJ) is at least ten percent of that due to the 
external field interacting with the thermal photon gas (see Eq. 
©). Above the dashed line, the number of pairs created in a 
typical strong optical laser four-volume Qi = 6 x 10 -26 cm 3 s, 
is larger than unity. 



We mention at this point that since our analysis assumes 
a thermal equilibrium, as soon as the first pair is created, 
the conditions under which our result is valid, are altered. 
Therefore, Eq. (|4|) should be understood in the sense 
of a probability per unit volume per unit time for pair 
creation to ensue, rather than a rate per unit volume. 

Inserting the asymptotic rate Eq. §5§ into Eq. (01, 
making the substitutions r\ = ui/T and y = 1 — cos#, 
with cos# = fc • x/<jj and noticing that at small x, only 
large values of r\ contribute to the integral, and we arrive 
at: 
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The remaining integrations can be also performed, giving 
the final result : 
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This result is valid if \f5 <C 1 and it is interesting for 
several reasons. The rate shows a non-perturbative de- 
pendency on both the temperature and the external field, 
in both its pre-exponcnt and exponent and in a way quite 
different from the classic Schwinger result quoted in the 
introduction. In Fig. [2] we plot the logarithm of the 
rate Pfk(T, Xe)i calculated from Eq. ([6]), times the four- 
dimensional volume X 4 = 1/m 4 = 7.4 x 10~ 53 cm 3 s. In 
the region above the dashed line, the number of pairs cre- 
ated in a typical optical strong laser beam four-volume 
^; = tiVi, where n = 10 fs, Vi — tt x (0.8 ^m) 2 x ctj = 



6 x 10 -12 cm 3 (where we temporarily recover the speed of 
light, c), is larger than one. Upon inspection of Fig. [21 we 
notice that for high enough temperatures, even with E = 
0.025, an exponential number of pairs can be produced 
in a typical optical laser pulse. This should be compared 
to the case of zero temperature, for which the expected 
number of pairs is identically zero. The solid line in Fig. 
[5]shows the significance of the background process of pure 
thermal pair creation (see Eq. ([!])). which dominates for 
large T/m and small xe- By comparing the analytical 
expressions in Eqs. ([1]) and <j6j> , we obtain the condition 
that the pure thermal pair production process is negli- 
gible if a(T/m)exp(-2m/T) < S 1 ' 4 exp(-4/v / 3l). The 
simpler but over-conservative condition for the validity of 
our method, T/m <C \/S <C 1, is obtained by comparing 
the exponents in Eqs. ([1} and ([6|). 

The physical line of argument given above also has 
a formal counterpart that yields additional insight into 
our result Averaging the one- loop polarization tensor 
over a thermal photon bath corresponds exactly to a two- 
loop calculation involving a zero-temperature electron- 
positron loop and a thermalised photon radiative correc- 
tion. In the finite-temperature Matsubara formalism, if 
one makes the assumption that only the photon propa- 
gator need be thermalised (justified below), the rate is 
given by: 



Pm(T, Xe) = -j 
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where k n = (2irinT, k). By using an integral relation 
to transform the sum into a complex integral (see Eq. 
(4.11) in (1^), one can show that Eq. J?]) and the ther- 
mal average Eq. ((3]), coincide. The two- loop effective 
action has been calculated in the Matsubara formalism 
for general constant- field backgrounds in [l9|, 26 1 in the 



form of a triple propertime-integral representation also 
involving a sum over the windings of the photon around 
the Euclidean finite-T cylinder. Taking the crossed-field 
limit, the result <j6j) arises from a saddle point of the 
propcrtime integrals, which occurs under the condition 
that 7t = T 2 /m 2 S <C 1. The requirement of small jt is 
not only in line with the assumptions made above, but 
also shows that our pair production density is dominated 
by non-perturbative Schwinger-type pair production in 
contrast to a multi-photon regime for large 7t- 

Another lesson to be learned from the effective action 
representation is that pair production in crossed fields is 
triggered by thermal photons with high winding number 
n in Euclidean space. Contributions to pair production 
arise from winding numbers n > 2-\/3 7t- In the 7t <C 1 
limit this is satisfied for all n > 1. In the general case, it 
shows that pair production requires a minimum amount 
of derealization of the virtual electron-positron pairs in 
order to acquire sufficient energy to become real. The 
winding number n is a measure for this derealization. 

Finally, we can also extract information about the fun- 
damental validity limit of the above reasoning. For in- 
creasing temperature, also the thermal fluctuations of 
the electron-positrons need to be taken into account. 
For T <C m, they are t ypi cally exponentially suppressed 
~ exp(- 



-m/T) 
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32j, such that we expect our ap- 
proximation of thermalising only the photon to hold for 
T< m. 

Our results also confirm the tendency observed for 
Schwinger pair production that the thermal contribution 
exceeds the vacuum contribution in the limit of weak 
fields j26|. In the present crossed-field case, this is partic- 
ularly evident, as the vacuum production rate in crossed 
fields is exactly zero. 

In order to justify the validity of using the asymptotic 
limit for single-photon pair creation, we can numerically 
calculate the ensemble pair density by inserting the full 
single-photon rate Eq. ([2]) into the thermal integral Eq. 
Q and compare it with the analytical expression in Eq. 
©. Upon doing so, we see that up to values of 8 = 0.03 
(the highest value of 5 plotted in Fig. [5] is 0.01), there is 
a maximum deviation of around 10%, as plotted in Fig. 
[3l Therefore, the asymptotic expression for the created 
pair density in Eq. <j6j) gives a good approximation to the 
exact value calculated by this method. We note that by 
passing to the integration variable 77 = uj/T in Eq. Q, it 
can be shown that the rate pth has the form pth(T, \e) = 
(T/m) 2 p t h{S) such that the explicit dependence on T/m 
cancels in \Ap t h(T, XE)\/ptb(T, xe)- 

We conclude by recapitulating the main results of the 
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FIG. 3. Relative difference \Ap t h{T, xe)\/ pth{T,XE) between 
the exact and asymptotic formulae for the thermal density of 
created pairs (using Eq. ([2]) in Eq. Q and using Eq. (B), 
respectively), plotted as a function of 5. 



letter. We have shown how, in a region of the thermal 
vacuum in equilibrium with an external constant crossed 
field, pair creation can indeed take place. By exten- 
sion, this conclusion also applies to regions of the vac- 
uum polarised by single plane waves, when in the tun- 
neling regime. The probability for this new process is 
non-perturbative in both the external field strength and 
the temperature of the thermal vacuum and is distinct 
from other analytical results for pair creation by external 
fields. For realistic field intensities orders of magnitude 
less than critical, a thermal background can seed the cre- 
ation of an exponential number of particles. Thus the old 
adage that the vacuum cannot be polarised by a single 
plane wave, although useful in understanding the the- 
oretical basis of vacuum polarisation, is an idealisation 
which yields when heat and hence any physically realis- 
tic scenario, is taken into account. 
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